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Abstract —This paper deals with the synthesis and design of 
impulse response digital systems exhibiting negative group 
delay (NGD). It is demonstrated that under certain 
condition, certain first order FIR and HR systems can 
provide the NGD function. Synthesis methods of NGD 
FIR/IIR transfer functions are introduced. Then, the basic 
properties specific to the NGD digital functions are 
established. To validate the theory proposed, NGD FIR and 
IIR systems were synthesized and implemented with Matlab. 
As expected, group delay frequency responses with NGD in 
baseband frequency well-correlated to the theoretic 
hypothesis were found. Then, thanks to the NGD effect, 
arbitrary form discrete transient signals with limited 
bandwidth propagating in time advance were demonstrated 
in time-domain. It is proved that this counterintuitive 
function does not forbid the causality principle and the 
NGD numerical system respects the functioning condition of 
the classical digital system. The NGD principle presented is 
potentially useful for the group delay correction processes 
and the signal pure delay reduction. 

Index Terns —digital system, negative group delay (NGD), 
synthesis method, NGD system design, delay equalizer 

I. Introduction 

Since the middle 1990s, theoretical and experimental 
investigations have demonstrated that certain electronic 
circuits are capable to operate with negative group delay 
(NGD) [ 1 ]-[4]. Similar to the classical electronic 
functions (filter, amplification, phase shifter...), it was 
found that the NGD circuits respect the criteria for the 
correct functioning of electronic systems. Under specific 
conditions between the NGD transfer function parameters 
and considered electrical signal bandwidth, it was shown 
several times that through the NGD circuit, the electrical 
signal output can propagate in time-advance of its input 
[5]-[9]. In this case, it was pointed out that this 
counterintuitive phenomenon does not forbid the 
causality principle [3], [4], 

At the beginning, the NGD electronic circuit was 
designed with low frequency components as the 
operational amplifier associated with passive networks 
[l]-[4]. Therefore, its applications were limited to very 
low frequency systems only up to some MHz. Then, 
NGD passive devices operating at microwave 
wavelengths were proposed [10]-[14], But analytical 
theories and experimentations confirmed that with such 
passive devices, the apparition of NGD effect is 
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systematically accompanied with significant losses [1], 
[10]-[13], 

To overcome this limitation, NGD microwave active 
circuit topologies were developed recently [15]-[17], It 
was stated that these NGD active topologies were 
susceptible to operate up to several GHz, offering then 
possibilities of loss compensations. So, several types of 
NGD circuit applications, in particular, in 
telecommunication area were introduced [15]-[17], 
Currently, one of technological issues which bottleneck 
the progress of current digital electronic systems is 
caused by the delay effects. As examples, it was reported 
that designers of microwave devices [18], [19] and 
electronic boards and circuits [20], [21] must take into 
account the different kinds of delays induced by the 
electrical links and interconnect buses. To cope with such 
aspects, various applications of NGD circuits based on 
the delay and phase compensation techniques were 
introduced recently. First, in the microwave area, low 
group delay broadband phase shifter [22] operating 
independently with the frequency were synthesized and 
implemented. In addition, enhanced structure of feed¬ 
forward amplifier was proposed [23]. Then, further 
application based on the group delay equalization 
dedicated to the analogue-numerical or mixed 
microelectronic systems was also established [24]-[27], It 
consists in the recovery of the mixed signal integrities 
degraded by the interconnect circuitries between the logic 
gates or electronic chips with possibilities of rise-/fall- 
time and propagation delay reductions [25]. 

In this paper, innovative principle of the digital 
systems exhibiting the outstanding NGD function is 
developed. Similarly to the aforementioned analogue 
delay compensation approaches, the NGD effects can be 
used for the diminution of the constraints caused by the 
fractional or pure delays in the digital filters and control 
systems [28]-[35], In addition, by transposing in 
numerical area the techniques introduced in [24], [25], 
the NGD digital systems could be employed for the 
equalization of the group delay based on the digital signal 
processing (DSP). 

For the better comprehension, this paper is structured 
in three main sections. Section II focus is on the analysis 
and synthesis of IIR and FIR digital systems capable to 
exhibit the NGD effect. Then, specific properties and 
characteristics of the NGD function will be also 
established. Section III is depicted with the verification 
and validation of the analytic concept established. Design 
examples and numerical experimentations run with 
Matlab illustrating the relevance of the theory will be 
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analyzed and discussed. Then, frequency-domain and 
transient analyses with different types of test signals will 
be made. Then, the last section is the conclusion of the 
paper. 


Synthesis Method'. At very low frequencies ( co « 0 ), 
the magnitude of H(jco) becomes: 

//(0) = • |(1+ £>)/(! +a)| (6) 


II. Theoretic Approach on the Elementary NGD 
Digital Impulse Responses 

First and foremost, let us consider the first order HR 
system excited by the discrete input signal and generating 
samples of output respectively, denoted x[n] and y[n] (n- 
th samples) at the discretized time t n =nT s (T s is the 
considered sampling time) defined by the difference 
equation: 


The same, the group delay frequency responses will be 
expressed as: 

T(0 ) = T s (b-a)/((\+a)(\+b)) (7) 

One finds that this quantity is negative under the 
condition: 

-—— < 0 < t = t > b — ci <0, (1 + 17)(1 + b) > 0 (8) 

(1 + £j)(l -b ^0 


y\n +1]= H 0 { x\n +1]+ b- x[«] }— a ■ y[«] (1) 

with // ( )>0, a and b (afb) are real constant. The transfer 
function canonical form of this system is written as: 

H(z) = Y(z)/X(z) = H 0 (l + b- z ~ l ) /(I + a ■ z _1 ) (2) 

Fig. 1 represents the block diagram of this HR system 
under consideration. 



To preserve the stability of the system understudy, 
along this paper. Id is supposed lower than 1. 


A. Analysis of Elementary NGD HR Filter 

It is well-known that the harmonic responses of the 
system expressed in Eq. (2) is realized replacing z by e ]c " 
( a>= 2/rf • T s is the angular frequency and j = V-T). So, 
the magnitude and phase responses are respectively 
expressed as: 


H(co) = If 


^/[l + b cos(ru)] 2 + b° sin 2 (®) 
-J[l +acos(®)] 2 +a 2 sin 2 (®) 


<p(co) = arctan 


b 2 sin° (co) 
1 + bcos(co) 


- arctan 


a 2 sin (co) 
1 + acos(co) 


(3) 


(4) 


Therefore, the group delay frequency response defined 
by r(co) = -d<p(co)/dco , is given by: 


t(co) = 


b[b + cos(fti)] a\a + cos(fti)] 


1 + 2bcos(co) +b 2 1 + 2acos(co) +a 2 


T s (5) 


In this case, an innovative type of HR system with 
negative group delay (NGD) is realized. By exploiting the 
previous frequency analysis, a synthesis method of NGD 
HR system can be established according the desired 
values of negative group delay r 0 and gain a at very 
low frequencies (co~ 0). Based-on the canonical form 
expressed in Eq. (2), the design of the NGD HR system 
depends on the determination of parameters H 0 , and b by 
fixing a (or inversely, b by fixing a). For example, by 
fixing the parameter a, the term b can be synthesized by 
inverting equation r(0) = r 0 . Therefore, the following 
synthesis formula is yielded: 

b = a ' T s + (1 + a ) T o (9 j 

T s -(1 + a)r Q 

It is interesting to note that in order to design an NGD 
HR system, r 0 must be negative. The parameter H 0 can 
be determined similarly, by solving the equation 
H( 0) = a . So, one gets the following expression: 


_ a(T s -(l + fl)r 0 ) 

(l + a)7; 


Knowing these synthesis formulae, it can be wondered 
what are the specific characteristics and properties of the 
HR systems with the NGD function. 

Properties of NGD HR Systems: In practice, the group 
delay of first order linear systems can be negative but its 
frequency bandwidth must be limited. The NGD cut-off 
frequencies can be determined from equation 
T{2.7tf c ■ T s ) = 0 . By using Eq. (5), it results the following 
formula: 


f c (k) = 


arccos(-(o +b)/(l + a -b))+ 2k7t 

2 tr s 


(ii) 


where A: is a constant integer. The targeted application of 
the NGD concept established in this paper concerns the 
field of base band signal processing. So, the operating 
signal frequency spectrum must start from zero to certain 
frequency limit (in order of GHz). For this reason, the 
first cut-off frequency can be assumed as the lowest of 
frequencies defined in Eq. (11), f=f(0). Moreover, 
substituting Eq. (9) into Eq. (11), the following relation 
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between the desired value of NGD r 0 and the cut-off 
angular frequency is extracted: 


f. =-arccos 

c 2 ttT s 


2a T s +(l-fl 2 )r 0 
(l-a 2 )r 0 -(l + o")^ 


( 12 ) 


This expression states that for the first order NGD HR 
systems, the cut-off angular frequency f c is inversely 
proportional to lr 0 l. In order to achieve a uniform 
propagation of each spectral component of the signal 
exciting the NGD systems expressed in Eq. (2), the cut¬ 
off frequency f. must be higher or equal to the signal 
bandwidth for the simplification here denoted /'. During 
the synthesis process, the maximal value of the desired 
NGD must be equal to: 

2a + (1 + a 2 ) cos(2cf s T s ) 

T °"“ ' (fl 2 -l)[l-cos(2^r f )] 

It is noteworthy that the performance of NGD systems 
depend on the surface A(r < 0) delimited by the group 

delay t(co) , x- and y-axis as will be seen in the next 
section about the numerical experimentation. This surface 
is mathematically defined as: 

co=2tf c T s 

Air < 0) = I" T(co)dco = <p( 0) - <p(2jf c ■ T s ) 
a >=o 


=> A{t < 0) = arctan 


(a-b)/yl(a 2 -l)(b 2 -l) 


(14) 


B. Analysis of Elementary NGD FIR Filter 

One can realize a FIR system from the HR transfer 
function expressed in Eq. (2) by taking a= 0. The 
canonical form becomes: 


H(z) = H 0 (l + b-z~ 1 ) 


(15) 


Frequency Analysis: The gain and phase harmonic 
responses are respectively expressed as: 

H(co) = H q --\/[l + Z?-cos(ft>)]“ +b 2 sin 2 (ft>) (16) 


<p(co ) = arctan 


b 2 sin° ( co) 
l+b- cos (eb) 


(17) 


Through the derivative calculation with respect to co, 
yields the group delay of the FIR system understudy: 


t{co) = 


b[b + cos(<y)] 


1 +2bcos(a>) + b~ 


(18) 


Synthesis Method: One can see that at very low 
frequencies, the magnitude and the group delay responses 
are written as respectively: 


H(0)=H () -\\+b\ 
r(0) = Th!(\ + b) 


(19) 

( 20 ) 


By considering \b\<\, this last expression can be 
negative under the condition: 

b/(l + b)<0 <=> (£> < 0,1 + £> > 0) (21) 

Similar to the previous case of the HR systems, the 
NGD cut-off angular frequencies can be determined with 
equation: 

cos[2 7tf c kT s ] = -b => f c = mCC ° S( _ b) (22) 

2nT s 

By inverting equation Eq. (20) and taking r(0) = r 0 , 
one finds the following synthesis relation: 

b ~ Tq l(T s — Tq ) (23) 

As r 0 must be chosen as to be negative, the absolute 
value \b\ must be lower than 1. Similarly, by solving 
equation H(0) = H , the second parameter of the NGD 
FIR proposed is established: 

Hq =H(1-t 0 ) (24) 


C. Analysis of High Order NGD Digital Function 

When cascading two first-order NGD FIR H\{jco ) and 
H 2 (Jco), a second order FIR system with the following 
transmittance is realized: 

H(jo)) = Hfjo))-H 2 (j of (25) 

So that, the total group delay is equal to: 

r(ty) = r, (co) + r 2 (co) (26) 

Meanwhile, the group delay value is proportional to 
the number of cells in cascade. Contrarily to the classical 
systems with positive delay, in the present case NGD 
systems, the negative delay increases with the number of 
cells. In time-domain, it corresponds to the increase of the 
signal time-advance instead of signal delay. This finding 
is important for the optimization of high order systems 
with significant value of NGD surface area in wide 
frequency band. By cascading k identical NGD digital 
systems, the total group delay must be equal to 
r k (co) = k ■ t(co) and the surface delimited by expressed 
in Eq. (14) the NGD curve must be 
A it (r<0) = k-A(r<0). 

III. Validation with Numerical Experiments 

The effectiveness of the previous theoretic concepts is 
intended to be verified in this section. Design examples 
of NGD HR and FIR systems are realized and tested via 
frequency- and time-domain numerical analyses. 
Different types of transient digital signals are considered. 
It is noteworthy that the numerical results presented in 
this section were computed with Matlab. 

A. Design and Synthesis of NGD FIR System 

As application of the previous concept, the desired 
values of magnitude and group delay are arbitrarily 
chosen respectively, H odB =-0.57 dB and r 0 = —2 ns. By 
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applying synthesis relations Eq. (23) and Eq. (24), the 
NGD FIR filter having the following difference equation 
is designed: 

>■(« +1) = 14.976{ x[n +1 ] - 0.937x[n] } (27) 

By implementing this NGD FIR function into Matlab 
program, the group delay plotted in Fig. 2 was obtained. 
As expected in theory, baseband NGD with minimal 
value equal to -2ns with cut-off frequency of about 
424MHz for the sampling time 7 V =0.13ns. This NGD 
frequency response is delimited by the shaded surface 
area shown in Fig. 2. So, an excellent agreement between 
the theory and the numerical experimentation is realized 
To illustrate the functioning of this NGD FIR, normalized 
discrete Gaussian input signals denoted X| and x 2 sampled 
in the time interval [0, 7’„ lat =40ns] with 301 sampling 
points with half-width at half maximal heights A7, =5.89 
ns and AT 2 = 1.47 ns are considered. 




Figure 2. Frequency responses of FIR expressed in (26): group delay 
(in top) andjft(xi t 2 ) (in bottom). 




Figure 3. Time-domain responses of NGD FIR expressed in (26) 
corresponding to inputs xi (in top) and X 2 (in bottom). 


As displayed in the bottom of Fig. 2, the frequency 
spectra of xi and x 2 present frequency bandwidths 
respectively of about f sl ~ 429MHz and / j2 ~l .71 GHz. 
Therefore, as can be seen in Fig. 3, the output signal y x 
practically similar to X| in time-advance of about -2 ns is 
obtained. However, as f s2 is widely higher than the cut off 
frequency of the NGD FIR understudy, the output signal 
y 2 (see Fig. 3) corresponding to the input x 2 is completely 
distorted. It is worth underlining that this time advance 
does not forbid the causality principle. As illustrated in 
Fig. 4, at the starting point of the signal, the output yi[0] 
appear at the same time as the input xi[0]. Then, as 
sketched in Fig. 4, thanks to the NGD function, this 
output behaves as an anticipated form of the input. 



10 15 20 

Time (nsec) 


Figure 4. Zoom in of transient responses of NGD FIR defined in (27): 
in top [0ns, 8ns] and in bottom [10ns, 20ns], 


Furthermore, as explained by the experimentation of 
truncated signal plotted in Fig. 5, the output value yi[kTJ 
is generated from X|[k7’ I ], xi[(k-l)7’J but not from the 
future value x 1 [k7’ s +lr 0 l]. Meanwhile, the future value of 
the input does not affect the present value of the output. 

2 


0 

7c 

g> "2 
CO 


-4 


-6 

0 10 20 30 

Time (nsec) 

Figure 5. Test of causality with truncated signal. 












- Input 
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The NGD FIR does not anticipate the signal 
discontinuity because this singularity point contains 
harmonic components outside the NGD bandwidth. It is 
important to note that the NGD function can operate also 
with different forms of smoothed input signals. As 
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demonstrated in Fig. 6, once again, output signals in 
advance compared to the input is occurred when injecting 
a signal with arbitrary form. 



Figure 6. Demonstration of time-advance with arbitrary signals. 


B. Design ofNGD HR System 

Similar numerical test as in previous subsection was 
made with an NGD HR system but in this case the desired 
values of gain and group delay are arbitrarily chosen 
equal to H 0dB = OdB and r () =—1.5 ns. By considering a 
pulse signal and 30ns periodical signal sampled in the 
time interval [0, r„ MV =60ns] with 501 sampling points. By 
applying synthesis relations Eq. (9) and Eq. (10), the 
NGD FIR system with difference equation was designed: 

y[n + 1] = 13.775{.r[n + lj- Q.921Ax[n ]}- 0.25y[«] (28) 



Figure 7. Time-domain responses of NGD HR filter expressed in (27). 

As explained in Fig. 7, a group delay response with 
NGD bandwidth of about 676MHz is realized with 
optimal value -1.5ns at very low frequencies. With 
Matlab numerical experimentations, once again, output 
signals in time-advance were confirmed by considering 
different types of inputs with arbitrary form as depicted in 
Fig. 8(a) and periodical form as displayed in Fig. 8(b). As 
aforementioned earlier, one of potential applications of 
the NGD digital systems lies on the principle of the delay 
compensation. Delay equalization process diagram is 
proposed in Fig. 9(a) by cascading a fractional delay 


having equation (rz 0 =int(lrol/73s-)) with the NGD HR 
defined in Eq. (27) as experimentation of the proposed 
technique efficiency. The test fractional delay was set at 
Time-Delay=\T () \= 1,5ns for verifying the delay 
compensation. As illustrated in Fig. 9(b), the delayed 
signal is given by x 0 (n.T s )=x(n.T s -DeIay). After 
interaction with the NGD HR system designed, as proved 
in Fig. 9(b), the NGD system output is approximately 
equal to the general input, y[n]~.r[u] for 
n>n 0 =mt(Delay/T s ). 



0 10 20 30 40 50 60 


Time (nsec) 



0 10 20 30 40 50 60 


Time (nsec) 

(b) 

Figure 8. Time-domain responses of NGD IIR filter expressed in (27): 
(a) arbitrary pulse signal - (b) 30ns periodical signal. 


4n\- 


? n ° 


x ( ln]=x{n-n ( } 


H h . 


z) ->y[n] 


(a) 



(b) 

Figure 9. (a) Delay equalization process diagram with NGD IIR 
system - (b) demonstration of transient signal compensation. 


IV. Concluding Remark and Future Work 

Several investigations on the NGD analogue devices 
based on the electronic analogue circuit were published in 
the literature [1]-[17], For this reason, innovative 
synthesis and design of IIR and FIR digital systems 
susceptible to exhibit NGD effect were devoted in this 
paper. By considering simple first-order transfer function 
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digital systems, it was demonstrated that NGD IIR/FIR 
functions can be realized. Then, synthesis method based 
on the desired value of NGD was established. Then, 
properties and characteristics specific to the NGD digital 
systems notably the determination of the NGD bandwidth 
and the effects of the association of NGD digital cells in 
cascade were developed. 

In order to substantiate the relevance of the concept 
established, NGD HR and FIR systems are designed and 
implemented into Matlab program. So, as expected, 
group delay response presenting NGD in base band 
frequency were obtained. Then, by considering different 
types of pulse signals, it was proved that when the 
frequency band of the input signal belongs to the NGD 
bandwidth, the output signal can propagate in time- 
advance to the input. It was confirmed also that the NGD 
function does not forbid the causality and the stability. In 
addition, it was demonstrated that this NGD effects can 
be potentially employed to the equalization of pure delay 
especially, in the numerical or mixed control or/and 
telecommunication systems. 

Due to the incessant increasing number of applications 
involving digital signal processing (DSP) and numerical 
filtering, the variety of requirements that have to be met 
by digital systems has increased considerably [36]-[39], 
Digital filters are integral parts of many DSP systems, 
including control systems, systems for audio and video 
processing, communication systems and systems for 
medical applications. In these systems, group delay 
equalizers play an important role by connecting in 
cascade with digital filters. In the continuation of this 
work, by improving the equalization method as described 
in [40]-[42], the study of the present NGD concept 
integration in the DSP for the synchronization of 
numerical data and delay correction in the complex high¬ 
speed system structures is in progress. 
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